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Abstract 

In this work we study the CP violation for J/ipKs and —> dAg up to 

leading and next-to-leading order QCD corrections in the standard model, two-Higgs- 
doublet model and the minimal supersymmetric extension of the standard model. We 
also study the effect of new physics on the branching ratios of these two decay modes. 
We find that within the parameter space constrained by the observation of the decay 
b S7, new physics does not affect the CP asymmetries greatly, and the prediction 
of new physics to the branching ratios of B^ J/ipKg and B^ —> ^Ag is the same as 
that of the standard model up to a minor discrepancy as far as the Yukama couplying 
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constants are perturbative. 
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I.INTRODUCTION 


The decay of J/ij^Ks is expected to be one of the most promising channels 

to stndy the CP violation in the B decays at the B factory. In the standard model CP 
violation is produced via the complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. Because the CKM matrix is unitary there is the unitarity equation 

= o {j^k). (1) 

i 

This equation can be represented as a closed triangle in the complex plane, which is 
called the unitarity triangle. Since the most poorly known entries in the CKM matrix are 
contained in the triangle equation 

VudV:, + VM + VuVl = 0 , ( 2 ) 

so this equation is most useful from the phenomenological point of view [P . This equation 
can be visualized as the unitarity triangle in Fig.l. The angles of a, /3, and 7 in Fig.l can 
be determined through the detection of the CP violation of the decay modes B —> tttt, 
B —> i^Ks, B 4>Ks, and Bg —> p^Kg respectively |]|, ||. In the standard model CP 
violation is only produced through the CKM mechanism. It is believed that if the sum of 
the three angles a, (3, and 7 is measured to be vr, the standard model gives the correct CP 
violating mechanism. If not, then there must be new physics beyond the standard model. 

The decay B^ —> J/'ll)Kg and B^ —> (l)Kg are based on the transition b qqs in the 
quack level, where q can be c and s. Their CP conjugate processes are B^ —> J/il^Kg 
and B^ —> (l)Kg respectively. Thus both the state and CP conjugate state of the initial B 
meson can decay into the same final state. It’s CP violating parameter, CP asymmetry, 
can be described byp 

^ /) - ^ f)]dt 

/o“[f(^pV(^) ^ /) + ns^git) - f)]dt 
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1 _ |g|2 _ 2Ima^mB/TB) , . 

(l + |C|2)[l + (AmB/rB)2]’ ^ ^ 

where ^ here ‘+’ for CP even |/), ' for CP odd |/) Q. For the present two 

decay modes that we study, ^/ = — (p)^ (p)ii: parameters p, q are dehned in the 

physical state of B and K mesons. The Kg meson is dehned as \Ks) = pk\K^) + Qk\K^), 
and the physical state of the neutral B meson is dehned as 

\Bl) = +(zb|bO), 

\Bh) = PB\B^)-qB\B^), (4) 

where L, H denote light and heavy respectively. In the standard model the decay B^ —> 
J/'ijjKg is dominated by the tree level contribution, and under the tree 

level approximation. The B^ —> <j)Ks is forbidden at the tree level, its relevant value of 
^ is also equal to —approximately. However, in some non-standard models such as 
the two-Higgs-doublet model (THDM) and the minimal supersymmetric extension of the 
standard model, CP violation can either be produced through the CKM mechanism, or 
through the Higgs sector, or from both. So in order to investigate the very nature of CP 
violation, the B decays should be studied carefully in both standard and non-standard 
models. 

In this work we study what impact the new mechanisms of non-standard model can 
give to the measurement of the CP violation of B^ —> J/'ijjKs and B^ —> ^Kg. To take 
into account the nonfactorization corrections, the inverse of the color number should 
be replaced with ;^ + X) where x stands for the nonfactorization corrections, and this 
replacement amounts to changing the color number from Nc = 3 H]. Here we consider 
three conditions with Nc = 2, 3 and oo. We also investigate the new physics effects on the 
branching ratios of B^ J/'ijjKg and B^ —> (l)Kg. The paper is organized as the following. 
In section II we study the CP violation and the branching ratio of B^ —> J/ipKg in the 
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SM, THDM amd the minimal supersymmetric extension of SM (MSSM). Section III is 
devoted to the study of —> cjiRg in those three models. Section IV is for the discussion 
and conclusion. 


II.The study of 

1) CP violation of J/'ijjKs in the standard model 

The low energy effective Hamiltonian relevant to our study is|^, 


'^eff 


Gf 

V2 


^ vJ Q?c'i(/x) + glc'2(/x) 


+ ^ QkCkilJ-) 

k=3 


+ n.c., 


(5) 


where Ck{^J-) (k=l, •••, 10) are Wilson Coefficients (WC) which are calculated in the 
renormalization group improved perturbation theory and include leading and next-to- 
leading order QCD corrections. Vq is the product of Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements and defined as 

{ v*^Vqb b ^ d transitions, 

VqsVqb b -I- s transitions. 

The ten operators are taken in the following form|^]: 



= {Uab(3)v-A{qpUa)v-A, 

Q 2 = {ub)v-A{qu)v-A, 

Qs 

= iqb)v-A'^iq'q')v-A, 

q' 

Qa = iqi3ba)v-A'^{q' aq'(3)v-A, 

q! 

Qs 

= {qb)v-AY^{q'q')v+A, 

9.' 

Q6 = {q/3ba)v-A'^iq'aq'p)v+A, (6) 

9^ 

Qi 

= \{qb)v-A^eq'{q'q')v+A, 

9' 

3 

Qs = i:{qi3bo)v-A^eqi{q' ^q'0)v+A, 

q' 

Qg 

= -j^{qb)v-A^eqi{q'q)v-A: 

3 

Qio = 2^q0K)v-A'^eqi{q'aq'/3)v-A, 
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where Qf and current-current operators, q=u, c. For q = c case, Qf and Q 2 

are obtained through making substitution u ^ c in and Qf- Qs ~ Qe ^^e QCD 
penguin operators, the sum ^ is runing over all flavors being active at fi = rrih scale, 
q' = {u, d, s, c, b}. Q-r ~ Qio are electroweak penguin operators, eq> are the electric charges 
of the relevant quarks in unit e which is the charge of the proton. The subscripts a, {3 are 
SU{3)c color indices. {V ± A) referes to 7^(1 ± 75 ). 

For our calculation we use the renormalization scheme independent form of the Wilson 
Coefficients 



= Cu 

C'2 = 

C2, 


= Cs-Ps/N,, 

c'4 

= C4 -|- Fs, 

a 

= C5 -Ps/Nc, 

C’o 

= Ce + F;, 

c^ 

= Cj + Pe, 

C' : 

= Cs, 


= Cg -|- Pe, 

C[o 

= ClO, 


where is[@, |, |] 

C(^) = (i -Fffas(/r)/47r-kfJaem(Ai)/47r) • C(/r), (8) 


and 


P, = 


^ - G{mq,q,fi) 


^ iVc / 


C^\ 


^ - G{mq,q,iJ,) 


G{mq,q, iJ.) = —4 / x{l — x)dxln 

Jo 


— x(l — x)q‘^ 


10 fJ- 

here q=u,c. The renormalization scheme independent Wilson coefficients C*i(^) has been 
given in Ref. jl^. 


Cl = -0.313, C 2 = 1.150, C 3 = 0.017, 
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C4 = -0.037, Cs = 0.010, Ce = -0.046, 

Cj = -0.001 -a, Cs = 0.049 -a, Cq = -1.321 • a, (9) 

Cio = 0.267 • a. 


With the effective Hamiltonian given in eq.(5) we obtain 






+ X! '^9 

q=u,c 


/^l r^l Q /~<l 

c' -I- — -I- r’' -I- — -I- —p (c' -I- — -I- r’' -I- 


Nr 


Nr 


Nr 


ij/i;\{cc)v-AmK^\{sh)v-Am, 


? 0 \ 


( 10 ) 


After we get eq.(lO), we are now in the position to calculate i^pKsi 


( 11 ) 


and 


_ vlb^td v*csNd 
\P) B \P) K VtbVl^ VcsV*a 

1 — p — it] 

I — p + ip 




where p, p are the parameters in the Wolfenstein parameterization for the CKM matrix. 
Substitute the numerical values of C- into eq.(lO) and (11), we can get the numerical result 
of C'lpKs ■ For the case without loop correction, only C[ and contribute, 


~ -e 


-2i0 


1 — p — ip 


1 — p + ip 

For the case with loop corrections all C[ - ■ ■ C^q contribute, 

^-2i/3 ll8-73 -p + i7? 


^pKs = -e 


118.73 — p — ip 

1 — p — ip 118.73 — p + ip 
1 — p + ip 118.73 — p — ip' 


( 12 ) 


(13) 
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Experimental results have constrained p and p into a small area in the p-p planeQ, approx¬ 
imately, —0.2 < p < 0.3, and 0.2 < p < 0.4. From eq.(12) and (13) we can see that loop 
diagram correction to is only in the order of 0(10“^). So the unitarity triangle angle 
P can be determined through the measurement of the CP asymmetry in J/^l)Ks up 

to an approximation of the order 0(10“^). The eq.(13) is obtained with = 3, when 
Nc = 2, oo the loop corrections are also in the same order of 0(10“^). This conclusion 
has been obtained through an estimate that the penguin contributions are propotional to 
VtbV-^l ~ A^, VcbV*g ~ and VubV*^ ~ A^ (A 0.22), so up to very small corretions the 
penguin contributions have the same weak phase as the tree diagram contribution, then 


the penguin contribution affects the CP violation extremely small [|l|, 1T|. Our calculation 
using the QCD corrected Hamiltanian confirms this conclusion. 


2) CP violation of in the two-Higgs-donblet model 

In the two-Higgs-doublet model |l^ the diagrams contributing to B^ —> J/ipKs are 
shown in Fig.2. The relevant low energy effective Hamiltonian has the same form as 
eq.(5), but the Wilson coefficients at the scale p = Mw contain an extra contribution 
from the charged Higgs loop diagrams. We calculated the new contributions at the scale 
p = Mw at first, and then using the two-loop-renormalization group equation to evolve 
them down to the scale p = mb- The initial conditions in the scale p = Mw are 

Ci{Mw) = Cf^{Mw) + (14) 

where the SM contribution have already been given in Ref. Q. In this paper, 

we calculate the contributions of charged Higgs loop diagrams, the results are 

ct{Mw) = 0 , 
ct{Mw) = 0 , 
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where 


ctiMw) 


Cf (Mw) 
ctiMw) 
ctiMw) 


as{Mw) 

247r 


Fh± {x't 


t)'! 


Stt ^ ^ 


as{Mw) 

247r 

as{Mw) 


Fh± {x't 


t)'! 


Stt 




C?^iMw) 


rv 

— [EH±{xt) - DH±{Xt)] , 


ctiMw) 


= 0 , 




a 

hvr 


Eh±{x'^ - ^H^{x'^ + 


1 


2sin?9. 


-DH±{Xt) 


C^o"{Mw) 


= 0 , 


x^ = 


2 

mi 


Eh± (x) 

Fh± (x) 

Dh± (x) 


i ?2 
^ U 

18 


i ?2 

^ u 
6 


38x — 79x^ + 47x^ Ax — 6x^ + 3x^ , 

■ H--71- Inx 


6(1 — x)3 


(1 — xY 


16x — 29x^ + 7x^ 2x — 3x^ 

m-x)^ + TT^'” 


i ?2 

-—xt 

2 


+ 


[l — X (1 — x)' 


rlnx 


(15) 


(16) 


where Fu is the coupling constant between the charged Higgs and the up- and down-type 
quarks. We consider two distinct two-Higgs-doublet modes, model I and model II, which 
naturally avoid tree-level flavor changing neutral current (FCNC). In model I, one doublet 
( 02 ) gives masses to all fermions and the other doublet (0i) essentially decouples from the 
fermions, and here F^ = cotp. In model II, 02 gives mass to the up-type quarks, while 
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the other {(pi) gives the down-type quarks masses, and in this model = cot(3, where 
tan(5 = V 2 /V 1 , which is the ratio of the vacuum expectation values (VEV) of the two Higgs 
doublet. 

After presenting the initial values of the Wilson coefficients Ci{M\y) in the two-Higgs- 
doublet model, now we shall evolve them down to the scale /x = mb, so that to take into 
account the QCD corrections. The renormalization-group equation for C{fi) is given byQ, 


d n, , d 


r^(^W 2 ^ '-Tf 2 

C{—^,g ,a )=7 {g ,a)C{- 






g‘^,a) 


(17) 


where 7 is the anomalous dimensions. A. Buras et al. have calculated the 10 x 10 anoma¬ 
lous dimension matrix involving current-current operators, QCD penguin operators, and 
electroweak penguin operators up to next-to-leading order QCD corrections in the stan¬ 
dard model|^. Now in the two-Higgs-doublet model and the supersymmetric extension 
of the standard model, the anomalous dimension matrix 7 is just the same as that in the 
standard model. In this section and the next we will use the anomalous dimension matrix 
7 calculated by Buras et al. to evolve the Wilson coefficients from scale fi = Mw down to 
the /X = mb scale. 

The solution of the renormalization-group equation (19) is givenQ by 


C{ii) = U{fi,Mw,a)C{Mw), 


(18) 


where U{fi, Mw, a) is the evolution matrix from Mw down to ^ < Mw- By resolving the 
eq.(17), we can get the numerical result for U{mb, Mw,Oi), which is presented in Appendix 

A. 

The difference between the phases of mixing in the standard model and in 

the two-Higgs-doublet model is trivial fl^j. There are two types of box diagrams which 
contribute to the B^ — B^ mixing: the diagrams with one W and one charged Higgs propa¬ 
gating in it and the diagrams with two charged Higgses. The charged Higgs contributions 
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to the mixing involve the same CKM factors and almost have the same phase as the W- 


box contributions]^, |^. So we will take {q/p)B = which is the same as in the 

standard model. The case of the mixing is similar to the mixing. The 

contributions of charged Higgs bosons to the mixing present a negligible phase 

difference from the standard model. 

With the Wilson coefficients calculated in the two-Higgs-doublet model, and using 
eq.(3), (10) and (11), we can get the results of the CP asymmetry of B^ J/il^Kg. We 
now discuss our numerical results. 

In the parameter space constrained by the observed rate of 6 —> sy, the possibility that 
charged-Higgs exchange induces large CP-violating effects through the complex phase of 


v^jvx has been ruled out|^, ^|. In the continuing work we will drop the complex phase 
of the VEV, i.e., only take into account the cases that U 2 /U 1 is a real number. 

The observation of the decay 6 —> sy gives a correlated bound to cot/3 and the charged 
Higgs mass Mii± [p^] 

\cot(5\{imGeV/MH±) < 0.8. (19) 


We make our calculation within this constraint of cot[5 and Mfj±, and also keep cot/3 < 3, 
which is the region of cot/3 that keeps the Yukawa couplying perturbative|0. We calculate 
three cases with the number of color Nc = 2, 3 and 00 . The result is that within the 
parameter space constrained by the experiment, charged Higgs effect reduces the CP 
asymmetry. But the effect is too small to distinguish the two-Higgs-doublet model from 
the SM. 


3) CP violation of B^ in the MSSM 

In the minimal supersymmetric extension of the standard model (MSSM), except for 
the usual particle fields in the standard model, there are the relevant supersymmetric 
partners, and an extra Higgs doublet. They add the new contributions to the process 
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of 6 ^ CCS. They are the contributions from: i)the charged Higgs loop diagrams, ii)the 
up-type squarks and chargino loop diagrams, iii) the down-type squarks, gluino and neu- 
tralino loop diagrams, which are shown in Fig.3. For the basic structure of the minimal 


supersymmetric extension of the standard model, we refer the readers to the Refs. |18, 19|. 

Because the assumption is quite strong that all the parameters in the soft supersym¬ 
metric breaking terms, such as the masses and couplings of the MSSM scalars and fermions, 
are unified at the grand unification (GUT) scale n = Mqjjt, we relax all the constraints 
to the parameters of MSSM from this assumption and investigate the effects of the more 
general structure of MSSM. So the free parameters are the following: the masses of the 
left and right-handed up and down-type squarks, the mixing angles of the 

right and left-handed squarks, a the wino and charged Higgs masses and M^±; the 
superpotential // and cot(3. For simplicity the following features should be assumed: 

(i) Supersymmetric loop diagram contributions to B ^ Xg'j decay are mainly from the 
charged Higgs and chargino exchange. 

(ii) The first two generations of up and down-type squarks are almost degenerate. 


The above features are inferred from the observations noted in Ref. [19|. With these 
simplicities the remaining free parameters are left to be: a common mass for the first 
two generations of the left and right-handed squarks; the top squark masses and 

the top squark mixing angle af, the masses and M ^±; the /r and cotf3. 

The minimal supersymmetric extension of the standard model gives the new contri¬ 
butions to the effective Hamiltonian of eq.(5) by adding some new terms into the Wilson 
coefficients. We express the new terms as (/i). Thus 




SUSY/ 


( 20 ) 


As in the SM, the Wilson coefficients are calculated in the scale /i = Mw at first, then use 
the renormalization group equation to evolve them down to the scale n = mb- The values 
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of have been given in Ref. Q. We calculate the supersymmetric contributions 

to the Wilson coefficients, 

at // = Mw scale. The authors of Ref. |^] have calculated the contributions of 
gluino diagrams to the process b —> q'qq. We refer to their results for the gluino diagram 
contribution to in this work. The initial conditions of 

the evolution are listed in the following, 


= 0 , 


(jSUSY 


= 0 , 


fiSUSY _ _ 

j = d S b ^9 


[sL{xj,y) - ^OL{xj,y)^ - -^Z 


“s rySUSY 


(jSUSY 


= - E 


a 


2 T\/f2 


■ ^ ,92 'n^a Stt 

j=d^s^b 9 


( 21 ) 


f-iSUSY ^ 2q:^ 

j=d,s,b ^2 g 


(sR{xj,y) - ^OR{xj,y)^ - 


rySUSY 


where |^] 


CSUSY ^ ) + |izsc/sv 

g^s,l> 92 SVT 


2 4 

SL{xj,y) = SR{xj,y) = -g{xj,y) + -f{xj,y), 


7 2 14 

OL = Q9{xj,y) - -f{xj,y) + -F{xj) + -f{xj), 


( 22 ) 


and 


17 14 

Or = --g{xj,y) + -f{xj,y) + -F{xj) + -/(x^). 


/(x,y) = 

g{x,y) = 


1 


X y 1 1 1 

mx — ^- TT^lny -- + 


y - X I (x - 1)2 (y - 1)2 

if" 


X — 1 y — 1 J ’ 


X y , 1 1 

mx — 7- TTTrlny -- + 


X - y 1 (x - 1)2 (y - 1)2 


X — 1 y — 1 
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F{x) = 

fix) = 


1 


(1 — x)^ 

1 n 


r 3 , 15 n 21 9,0 o 1 

< “2®'^ ~2^ ~ ~2^ 2 ~ + 3x + l)lnx j-, 


(1 -x)4 i3 


Xj = mj/rfig, y = rhyrh 


The values of , (i = 7, • • •, 10) are 


! 3 2 

— —X + 3x 

—-— Inx 
0 

- 2 / ~ 2 

= 


, 10 ) are 


(jSUSY _ 

“ vSUSY 

Gvr ^ ’ 

riSUSY 

Ug — 

0 , 

r<SUSY 

09 — 

“ vSUSY 
7~-'2 ’ 

DVT 

r'SUSY 

*-^10 — 

0 . 


(23) 


(24) 


And in eq.(21) and (24) are presented in Appendix B. 

QCD corrections up to the next-to-leading order are taken into account by evolving the 
initial values of the Wilson coefficients down to the scale y = mb- Because the anomalous 
dimension matrices of the ten operators, C* (i = 1, • • •, 10), are the same as in the SM, 
we can still use the evolution matrix U{mb, M\y, a) listed in Appendix A. We perform the 
evolution using eq.(18) as in section 2 ). 

After obtaining the Wilson coefficients Ci {i = 1, • • •, 10) with the QCD corrections 
up to the next-to-leading order in MSSM, we can still use eq.(3), (10), (11) to calculate 
the CP violating parameter with the Wilson coefficients of the MSSM. We perform those 
calculations within the parameter space allowed by experimental limit. For example, the 
experiments give the lower bounds of squark and chargino mass are 176 GeV and 45 GeV, 
respectively )^] . 

Because the complex phase of mixing is the same as the standard model up 

to the minor correction of the order of (ruc/mt)^ or less[^], we do not pay much attention 
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to it here. 


We performed our calculation within the very large ranges of all the parameters (m^, 
H, rrig, rriuL, '^i) which are allowed by the experiment. We find that if we 

keep the Yukawa couplying is perturbative {cotj3 < 3)|l^, the discrepancy from the SM 
is extremely small, it only happens in the third desimal number after the zero point. 

4) The branching ratio of J/'ijjKs in the SM, THDM and MSSM 


Because there are great theoretical uncertainties in the calculation of the branching 
ratio of —> J/^|JK^ in this section we do not want to give an precise prediction of the 
branching ratio. We only want to give a comparison of the predictions of these three 
models, and to probe the possibility to distinguish them by measuring the branching ratio 


of this decay mode. We calculate the branching ratio according to the BSW method |^|. 
In our calculation we take = 0.386GeV, Ff^{0) = 0.379 and Ff^{q'^) = 
where = 5.43GeV. The result of our calculation with the color number being taken 
as Nc = 2, 3 and oo is that the difference between the predictions of the SM and MSSM 
is small, it is only up to a few percent level. 


Ill The Study of (t)Ks 


In this section we report the results of our study of B^ (j>Ks using the effective 
Hamiltonian given in eq.(5), (14), (18) and (20). Because there is no tree level diagram 
contribution to B^ —> in the SM, new physics may have an observable effect. Using 

the effective Hamiltonian given in the last section, the amplitude of B^ is 




Gf . V- 

V ^ q=u,c 


(1 + + (1 + + G'^ + —Cq + -Cs 


Nr 


Nr 


Gj + wC’s 
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+(1 + + l)C'io I {4>\iss)v-A\0){K^\isb)v-A\B^) 


(25) 


where (i = 1, • • •, 10) should be the renormalization scheme independent Wilson Coef- 
hcients calculated in the SM, THDM and MSSM, respectively, e* = — ^ is the charge of 
the s-quark. 

The Decay width of a meson at rest decaying into cj) and Kg is 


T{B^ ^ cj^Kg) = j-\{cj,Kg\Heff\B^)\^^, 

B 


where 


\P\ = 


{{Ml - (M^ + MkMmI - {M^ - 

‘IMq 


(26) 


(27) 


is the momentum of cj) or Kg. The corresponding branching ratio is given by 

T{B^ ^ <t>Kg) 


Bf,r{B^ ^(t>Kg) = 


^ tot 


(28) 


We take r£° = 4.22 x 10“^^GeV |21] in our calculation. The hadronic matrix ele¬ 
ment (i?i)|(ss)v'-A|0)(K‘^|(s?))y-A|.B°) is calculated in the Bauer, Stech, and Wirbel (BSW) 
method [^] . The one-body vector matrix element of (V-A) current is 


((/>|(ss)y_^| 0 ) = 


(29) 


where = 0.233 GeV, M^f, = 1.02GeV |^|, is the polarization of (j) meson. The 
two-body pseudoscalar-pseudoscalar matrix element of the vector current is 1 ^ 


{K^\is^^{l-^,)b\B^) = (^Pb + Pk - FBK{q^l-) + ^^^^^q>.FBKiq^0+), 

(30) 

2 

where q = Pb — Pksi FBK{q^-, 1 ~) = Fbk{0, 1~)/(1 “ m^”)- calculation, we take 


Fbk{0A ) = 0.379, Ml- = 5.43GeV [23|. Now with eq.(3), (25), and (28), we can 
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calculate the CP asymmetry and branching ratio of Our calculation shows 

that if we take cot(3 < 3 |17], the difference between the preictions of new physics (THDM, 
MSSM) and SM is too small to be detectable in the future B factory. 


V. Conclusion and discussion 

We studied the CP asymmetry and the branching ratio of —> J/ipKs and B^ —> 4)Ks up 
to the leading and next-to-leading order QCD corrections in the standard model, the two- 
Higgs-doublet model and the minimal supersymmetric extension of the standard model. In 
our calculation we can neglect the charged Higgs tree diagram contribution in the b ccs 
decay. Because this diagram gives the contribution to the operator C 7 ^(l -|- 75 ) 057^(1 — 


75 ) 6 , whose coefficient is —2cot^f3 jl^f , and eq.(19) means cot^/3 < ( 




O.Srric ^2 
. lOOGeyl • 


This 


u± " h ± 

implies the contribution of the charged Higgs tree diagram is too small, and can be safely 

neglected. So in the two-Higgs-doublet model and the MSSM, only some loop diagrams 
are added. Because of the coupling constant a and Og suppression, the W-boson tree 
diagram dominates. Thus there will be no large difference between the CP asymmetries, 
branching ratios of B^ —> J/%l)Ks predicted by these three models. 

For B^ —> (fiKg decay, new physics only contributes an overall strong phase, and it 
does not affect the weak phase much. So in the THDM amd MSSM, new physics has no 
large effect in the CP asymmetry of B^ —> (pKg. 

After finishing our work, we noticed that the next-to-leading QCD corrections to H —> 
Xg^ decay are calculated in the THDM in Ref. 




Finally we come to our conclusion: (i) The tree diagram calculation for Acp in B^ —> 
J/'ijjKg is reliable, which is independent of the nonfactorization corrections, (ii) It is 
difficult to distinguish the SM from the two Higgs-doublet model and MSSM by measuring 
CP asymmetries Acp and the branching fraction of B^ J/'ijjKg and B^ —> (jjKg at B 
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factories. 
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Appendix A 


The evolution matrix from Mw down to /x = mb scale, where we take mb = 5.0GeV. 

U{mb, Mw, a) = 


1 . 115 , 

- 0 . 246 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 

- 0 . 246 , 

1 . 12 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 

- 0 . 007 , 

0 . 012 , 

1 . 10 , 

- 0 . 203 , 

0 . 030 , 

0 . 078 , 

0 . 003 , 

0 . 008 , 

- 0 . 012 , 

0.012 

0 . 006 , 

- 0 . 033 , 

- 0 . 285 , 

0 . 977 , 

- 0 . 017 , 

- 0 . 169 , 

- 0 . 002 , 

- 0 . 017 , 

0 . 034 , 

- 0.018 

0 . 003 , 

0 . 009 , 

0 . 035 , 

0 . 043 , 

0 . 903 , 

0 . 096 , 

- 0 . 004 , 

0 . 004 , 

- 0 . 007 , 

0.001 

0 . 005 , 

- 0 . 038 , 

- 0 . 059 , 

- 0 . 172 , 

0 . 317 , 

1 . 717 , 

- 0 . 003 , 

- 0 . 029 , 

0 . 039 , 

- 0.024 

- 0 . 003 , 

- 0 . 001 , 

- 0 . 001 , 

0 . 003 , 

- 0 . 006 , 

- 0 . 001 , 

0 . 913 , 

0 . 053 , 

- 0 . 0096 , 

- 0.004 

- 0 . 0008 , 

0 . 0001 , 

- 0 . 0003 , 

0 . 001 , 

- 0 . 0026 , 

- 0 . 0104 , 

0 . 326 , 

1 . 949 , 

- 0 . 0028 , 

- 0.0005 

- 0 . 004 , 

- 0 . 0007 , 

0 . 0038 , 

0 . 0016 , 

- 0 . 0019 , 

- 0 . 001 , 

- 0 . 010 , 

- 0 . 004 , 

1 . 101 , 

- 0.2497 

0 . 0003 , 

- 0 . 0001 , 

- 0 . 001 , 

0 . 005 , 

0 . 0002 , 

- 0 . 0001 , 

0 . 001 , 

0 . 0009 , 

- 0 . 244 , 

1.112 


(dl) 
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Appendix B 


The functions and appearing in the Wilson coefficient Cj.^g^ 


in eq.(21) and (24). Our result is consistent with Ref. |1£] after replacing the lepton 
charge with the relevant quark charge. 

1) contribution from charged Higgs loops with a z-boson coupling to the quark pair; 

1 


Yi = -cor/3xtfi 


9 

mf 


m 


H±. 


^2 = - 2] COt^lSxtfl^ 


4 \siv?9,, 


m 


H±. 


(HI) 


2) contribution from charged Higgs loops with a photon coupling to the quark pair: 

1 


mr 


3) contribution from charged Higgs loops with a gluon coupling to the quark pair; 

1 


(H2) 


Z = -cot 13 h 


mi 

TI. 


(B3) 


4) contribution from chargino loops with a z-boson coupling to the quark pair: 


6 2 


Hi = - 


92'^ts'^tb j 




H2 = 


-C2(?T^L,m|±,m|±)5Ai?47iVji + -m-±m-±co(m|^,m|±,m|±)5Ai?C/nC^ji } , 


6 2 




siin?9. 






-C2{ml^,m^±,rn^±)5 abViiVji + -m~±m-±co{rnl^,mi±,mi±)6ABUiiU}^ } . 


m 
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5) contribution from chargino loops with a photon coupling to the quark pair: 


Vi = V2 = -- 


1 




2 

mt± 
Xi 

^Ia 


QgMsVtb 

6) contribution from chargino loops with a gluon coupling to the quark pair: 

1 


Z = -- 




2 

mt± 

Xi 


m. 


UA 


iB5) 


iB6) 


^aMsVtb 

7) contribution from neutralino loops with a z-boson coupling to the quark pair: 

-C2(m|^,m|o,m|o)(5AB(iV/*3iVj3 - N^^Nj^) - ^m^om^oco(m|^,m|o,m|o)(5 ab(A^/3A^j3 


-NuN*j^)}, 


Y2 = 


siv?6,. 


- 2 


6 4 


^ 2 „,* 


/ 92'^ts^tb 


E E (^E)i4(^E)B3{c2(m|o,m|^,m|j(r^«r^«+)AB<5/j 




-NmN},)}. 


m 


8) contribution from neutralino loops with a photon coupling to the quark pair: 

6 4 


Yi=Y2 = 


1 




2 

mlo 

Xi 


m- 


^^92V*tsVtb 

9) contribution from gluino loops with a z-boson coupling to the quark pair: 


(B8) 


h2 = 


8ffi 


^g^VtsUtb \sin‘^e. 


- 2 


E 

A,B=1 


(r^^)+^(r^^)s3C2(m?,m|^,m|j(r^«r^«+)AB. 
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{B9) 


10) contribution from gluino loops with a photon coupling to the quark pair: 


Yi=Y2 = 






m~ 


8lg2vtsVtb 

The functions appearing in the eq.(Bl)- • - (BIO) are given by 


m- 


dA , 


h{x) 

f2{x) 

hix) 

fiix) 

f5{x) 

hix) 

,2 2 2 


+ 


jinx, 


1 — X (1 — 

38x — 79x^ + 47x^ 4x — 6x^ + 

6(1 — x)3 (1 — x)‘^ 

16x — 29x^ + 7x^ 2x — 3x^ 


Inx, 


6(1 — x)3 (1 — x)^ 

52 — lOlx + 43x^ 6 — 9x + 2x^ 


+ 


+ 


(1-x) 


-Inx, 


C2{mi,ni2,ni^) = --- 


6(1 — x)3 

2 — 7x + llx^ 

6(1 — x)3 ' (1 — xY 

2 — 7x + llx^ 6x^ , 

-^-1- T- Tjlnx, 

(1 — x)3 (1 — x)^ 

/ m\ln{m\/ 

I (mf — rnYj{m\ — m'^) 

3 1 / m\ln{m1/fj?) 


Inx, 

^ Inx, 


(BIO) 


+ (mi ^ m 2 ) + (m-i ^ m3) , 


+ (mi <-> m2) + (mi <-> m3) 


8 4 I (mf — m2)(mi — m3) 

The symbol conventions in eq.(Bl- • • BIO) are the same as in Ref. [p^]. 
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Figure Captions 


Fig.l; The unitarity triangle of CKM matrix in the complex plane. 

Fig. 2; Diagrams contributing to the process b ccs in the two-Higgs-doublet model, 
which mediate J/^|JKs decay. 

Fig.3; Diagrams contributing to the process b ccs in the MSSM. 
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